Jet Charge at the LHC 
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Knowing the charge of the parton initiating a light-quark jet could be extremely useful both for 
testing aspects of the Standard Model and for characterizing potential beyond-the-Standard-Model 
signals. We show that despite the complications of hadronization and out-of-jet radiation such as 
pile-up, a weighted sum of the charges of a jet's constituents can be used at the LHC to distinguish 
among jets with different charges. Potential applications include measuring electroweak quantum 
numbers of hadronically decaying resonances or supersymmetric particles, as well as Standard Model 
tests, such as jet charge in dijet events or in hadronically-decaying W bosons in tt events. We 
develop a systematically improvable method to calculate moments of these charge distributions by 
combining multi-hadron fragmentation functions with perturbative jet functions and pertubative 
evolution equations. We show that the dependence on energy and jet size for the average and width 
of the jet charge can be calculated despite the large experimental uncertainty on fragmentation 
functions. These calculations can provide a validation tool for data independent of Monte-Carlo 
fragmentation models. 



The Large Hadron Collider (LHC) at CERN provides 
an opportunity to explore properties of the Standard 
Model in unprecedented detail and to search for physics 
beyond the Standard Model in previously unfathomable 
ways. The exquisite detectors at atlas and CMS let us 
go beyond treating jets simply as 4-momenta to treating 
them as objects with substructure and quantum num- 
bers. A traditional example is whether a jet was likely 
to have originated from a 6-parton. At the LHC, one 
can additionally explore whether a jet has subjet con- 
stituents, as from a boosted heavy object decay [l|, Q, 
or whether it originated from a quark or gluon Q. See 
Ref. 0] for a recent review of jet substructure. Here we 
consider the feasibility of measuring the electric charge 
of a jet. 

The idea of correlating a jet-based observable to the 
charge of the underlying hard parton has a long his- 
tory. In an effort to determine the extent to which jets 
from hadron collisions were similar to jets from leptonic 
collisions, Field and Feynman [f| argued that aggregate 
jet properties such as jet charge could be measured and 
compared. The subsequent measurement at Fermilab @ 
and CERN Q in charged-current deep-inelastic scatter- 
ing experiments showed clear up- and down-quark jet 
discrimination, confirming aspects of the parton model. 
Another important historical application was the light- 
quark forward-backward asymmetry in e + e~ collisions, a 
precision electroweak observable [8j. Despite its histori- 
cal importance, there seem to have been no attempts yet 
at measuring the charge of light-quark jets at the LHC. 

Most experimental studies of jet charge measured vari- 
ants of a momentum-weighted jet charge. We define the 



pr- weighted jet charge for a jet of flavor 



i as 



Q l = 



(1) 



where the sum is over all particles in the jet, Qj is the 
integer charge of the color-neutral object observed, pip 
is the magnitude of its transverse momentum and k is a 
free parameter. A common variant uses energy instead 
of pr- Values of k between 0.2 and 1 have been used in 
experimental studies d, [|| . 

In hadron-hadron collisions at high energy, such as at 
the LHC, the particle multiplicities in the final state are 
significantly larger than at low energy and at e + e _ or 
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1. Distributions of Q l K for various parton flavors with 
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Ncvcnt = 50 Dijet composition 




FIG. 2. Distinguishing W' from Z' with a log-likelihood dis- 
criminant, for different values of n. Even with only 50 events 
the samples are extremely well separated. 



lcpton-hadron colliders. Thus, one would expect that 
measuring the charge of a light-quark jet at the LHC 
should be difficult, with the primordial quark charge 
quickly getting washed out. However, this turns out not 
to be the case. For example, Fig. [T] shows distributions 
of Q l K for u, u, d, d and g jets for two values of k. One can 
clearly see that Q l K will be useful for measuring jet charge. 
Moreover, as we will show, the energy and and jet-size 
dependence of moments of jet-charge distributions can 
be calculated in perturbative QCD. 

To get an impression of how much data is needed for 
Q l K to be useful, we consider measurements designed to 
distinguish charged from neutral vector resonances. To 
be concrete, we consider scaled-up W and Z bosons at a 
mass of 1 TeV decaying into light quark jets Q. Simply 
cutting on the sum of the Q l K of the hardest two jets 
in each event we can distinguish the two samples with 
95% confidence using around 30 events. We find that 
the best discriminating power is achieved for k ~ 0.3. 
A more sophisticated log-likelihood discriminant based 
on the two-dimensional jet charge distribution is shown 
in Fig. [21 where ~ 4er separation of the two samples is 
achievable with 50 events. 

For another phcnomenologically relevant application 
of jet charge consider a simplified supcrsymmctric model 
with squarks pair produced through i-channcl gluino ex- 
change and decaying as q — > q + Xa^ At iriq = rrig = 
1.5 TeV such a model is still allowed [lOj, although it will 
come under scrutiny with the next round of 8 TeV data. 
Due to the high concentration of up-type valence quarks 
at large x, the di-squark production process yields many 
events with two hard up-type jets and missing energy, in 
contrast to the background (dominated by F+jets) where 
the two hardest jets are rarely both ups. Adopting a set 
of cuts similar to those of Rcf. Ildl, we estimate if an 



FIG. 3. Final state composition in dijet production. 




400 600 800 

Dijet mass [GeV] 



FIG. 4. Sum of the two jet charges in dijet events, for various 
k. The growth with dijet invariant mass reflects the larger 
fraction of valence quark PDFs at large x and corresponding 
decrease in gg final states. 

excess is seen in 2 jets and missing energy channel, the 
increased concentration of up quarks could be measured 
above the 2a level with 25 fb _1 of 8 TeV data, provid- 
ing unique insights into the flavor structure of the new 
physics. 

To trust a measurement of jet charge, it is important 
to test it on samples of known composition. While pro- 
ton collisions do not generally provide clean samples of 
pure up- or down-quark jets, there are still ways to vali- 
date the method on data. For example, dijet production 
has an enormous cross section at the LHC and the frac- 
tion of jets originating from different partons is directly 
determined by the parton distribution functions (PDFs) . 
At larger energies the valence quark PDFs dominate over 
gluon or sea quark PDFs, producing more charged final 
states, as can be seen in see Fig. [3] The mean total jet 
charge in dijet events is shown in Fig. U for various values 
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Distribution of jet charge for Ws from top pairs 
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FIG. 5. Sum of jet charges of the two non b-jets in semi- 
leptonic ti events with a positively (solid) or negatively 
(dashed) charged lepton. 

of K. Verifying the trend in this plot on LHC data would 
help validate jet charge. 

Another sample of interest for validating jet charge is 
hadronically decaying W bosons coming from top decays. 
In a semi-leptonic ti sample, the lcptonically decaying 
W can be used to determine the two charges of the jets 
from the hadronically decaying W . The distributions of 
these charges can then be compared to expectations, an 
example comparison is shown in Fig. [5] Validating this 
simulation on data would establish weighted jet charge 
as a trustworthy tool, which could then be used for new 
physics applications. Perhaps it could even be employed 
within the context of W decays to help with top-tagging 
or W polarization measurements. 

Next, we consider the effects of pile- up and contam- 
ination on jet charge. One might worry that at high 
luminosity jet charge would be diluted by pile-up events, 
as up to 0(100) proton-proton collisions can take place in 
the same bunch crossing. However, the products of these 
interactions tend to be soft, and are thus assigned little 
weight as long as k is not too small. Further, charged 
particles can be traced to their collision vertex allowing 
most contamination to be removed. Finally, jet grooming 
techniques like trimming (llj can be applied to further 
reduce contamination. Wc present a comparison of ef- 
fects of contamination and techniques to mitigate it in 
Fig. EI 

Having demonstrated the practicality of jet charge for 
new physics searches and proposed ways to validate it on 
standard model data, we now turn to the feasibility of 
systematically improvable jet charge calculations. While 
Monte-Carlo programs like pythia often provide an ex- 
cellent approximation to full quantum chromodynamics, 
they are only valid to leading-order in perturbation the- 
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FIG. 6. Comparison of W' vs. Z' discrimination subject 
to contamination from initial state radiation (ISR), multiple 
interactions (MI) , and pile-up events. We also show the result 
with and without jet trimming (i? su b = 0.2, / cut = 0.02). 

ory including the resummation of leading Sudakov double 
logarithms [121 ] . 

A precise calculation of jet charge is challenging be- 
cause it is not an infrared-safe quantity. Jet charge is 
sensitive to hadronization and cannot be calculated with- 
out knowledge of the fragmentation functions Dj(x,fi). 
These functions give the average probability that a 
hadron h will be produced by a parton j with the hadron 
having a fraction z of the parton's energy. Fragmentation 
functions, like parton distribution functions, are non- 
perturbative objects with perturbative evolution equa- 
tions which simplify in moment space. The Mellin mo- 
ments are defined by 

5>,/x)= / dxx v D^(x,n), (2) 
Jo 

which evolve through local renormalization group equa- 
tions, just like the moments of parton distribution func- 
tions. 

We first consider the average value of the jet charge 
«£> = — [d*Q? K = Uzz^Qh—^^, (3) 

Oj e t J ' h jot 

where z = Eh/E- iet is the fraction of the jet's energy the 
hadron carries. For narrow jets z ~ Pt/Pt ■ 

To connect to the fragmentation functions, we first ob- 
serve that for k > the the charge is dominated by 
collincar and not soft radiation. Thus the contributions 
of the hard and soft sectors of phase space, while con- 
tributing to the formation of the jet, should have a sup- 
pressed effect on Q* . We can therefore use the fragment- 
ing jet functions introduced in Refs. [HI, E3] to write 

1 dcrfrgjet 1 f 1 dx Jjj (E, R, f , ft) h 

<7jet dz 16tt3^W 2 x Ji(E,R,n) j{, ^ h 

(4) 
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FIG. 7. Comparison of theory prediction for the average jet 
charge (bands) to pythia (squares and circles for d and u 
quarks). Normalizing to 1 at E = 100 GeV and -R = 0.5 
removes the dependence on the nonperturbative input and 
quark flavor. 

Here Ji {E, i?, /z) is a jet function and Jij {E 1 R, x, /i) a set 
of calculable coefficients which depend on the jet defini- 
tion and flavor i of the hard parton originating the jet. 
The hard and soft contributions conveniently cancel in 
this ratio. Therefore 

with related to Jij by a Mellin-transform as in 
Eq. ©. By charge conjugation J2hQhE>g(n, fj.) = 0, 
so in particular (Q 9 K ) = 0. We have checked that the 
/j-dependence of Jij/Ji exactly compensates for the [i- 
dependance of the fragmentation functions at order a s . 

We have written both Ji(E,R 7 fj,) and Jij(E, fj.) 
as if they depend on the energy E and size R of the 
jet, however, these functions only give a valid description 
to leading power of a single scale corresponding to the 
transverse size of the jet. Here we use the e + e~ version 
of anti-fey jets of size R, for which the natural scale is 
fij = 2£'tan(i?/2) [l5j] . We can therefore calculate the 
average jet charge by evaluating the Mellin-moments of 
fragmentation functions at the scale fij and multiplying 
by the jet functions. 

Since only one linear combination of fragmentation 
functions appears in Eq.([5]), the theoretical prediction 
is not significantly limited by the large uncertainty on 
Dj(n,[j,). One can simply measure £)^(k, /i) by observ- 
ing the average jet charge for each flavor at one value for 
fi and then using the theoretical calculation to predict 
it at other values. In the absence of data, we simulate 
such a comparison using pythia. The result is shown in 
Figure [7] for various values of k and R, and normalized at 
a reference point. Already we can see a clear agreement 
between the theory and pythia. 

To calculate other properties of the jet charge dis- 
tribution requires correlations among hadrons. For ex- 
ample, we can consider the width of the jet charge, 



FIG. 8. Comparison of theory prediction (bands) for the 
width of the jet charge distribution to pythia (squares and 
circles for d and u). We take pythia at E — 100 GeV and 
R = 0.5 as input. 



(T J J 2 = {Q l K y - {{Q l K ) 2 ). This depends on the moment 
<(Qk) 2 ) = E E Jdz 1 ---dz n (Q 1 z? + ... + Q„z-) 2 

(6) 



n hi,...,h n ' 

1 d fl Q- fel .../ lTi gj ot 

o-jot dzi ■ ■ ■ dz n 



where the sum runs over all hadronic final states. After 
integrating over most of the Zi and including a factor of 
i for identical hadrons, this simplifies to 



•J I, 



1 da h 



Gjet 



Ojet dz 



[dz 1 dz 2 ZiZ$ V Qh^Qh 2 — 

hi,h 2 J 



(7) 

1 dah 1 h 2 ejct 



dzi dzi 



The first term on the right hand side can be expressed 
in terms of products of fragmentation functions and jet 
functions as for {Q l K )- The second term can be expressed 
in terms of something we call a dihadron fragmenting jet 
function, 2 . Its matching onto (dihadron) fragmen- 
tation functions is given by 



G^ h2 (E,R : z 1 ,z 2 ,fi) 



(8) 



hih 2 f z l z 2 



_■ 7, J 



The second term is due to a perturbative parton splitting 
before hadronization and only starts at 1-loop order, 

J$>(E,R,u,v,n) - j£\E,R,u,u)S(l-u-v)5 kMij) , 

(9) 
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where Sk, a (ij) indicates that the flavor k is completely 
fixed by ij. E.g. a(qq) = g, a(gq) = q. We then find 



{(Ql) 2 ) 



E 



16tt 3 ^ J q {E,R,n) 

hi ,h 2 



h 



(10) 



(For a gluon jet, which we do not consider here, there is 
a contribution from the last line of Eq. ([8]) corresponding 
to a perturbative g — > qq splitting.) We have checked 
that this equation is RG invariant at order a s . 

Unfortunately, the dihadron fragmentation functions 
are even more poorly known than the regular fragmenta- 
tion functions. However, we can use the same trick as for 
the average jet charge to calculate the E and R depen- 
dence of the width, given measurements at some reference 
scale. As with the average jet charge, we can now calcu- 
late the width by fitting one parameter for each flavor, 
corresponding to the term in brackets in Eq. (jlOp . and 
predicting the E and R dependence. Results compared 
to pythia for the width are shown in Fig. [5] and show 
good agreement. We did not include the gluon mixing 
contribution, which will increase the width and possibly 
reduce the agreement with pythia. 

To go beyond the average and the width, for example 
to the 3rd or higher moments, multi-hadron fragmenta- 
tion functions would be needed. From a practical point 
of view, such functions are nearly impossible to measure 
with any precision. However, we have found that dis- 
criminating power of jet charge is nearly as strong with 
just the average and width as it is with the full differen- 
tial jet charge distribution. It follows that accurate cal- 
culations of the phenomenologically relevant part of jet 
charge distributions are achievable with the formalism 
we have introduced in this paper. The full fragmenting 
jet functions, both for the single hadron and dihadron 
case, and the evolution kernels, are now known at 1-loop 
order. To see whether higher precision is required, and 
to explore the importance of power corrections, requires 
some LHC data to compare with. 

As we have shown, the weighted jet charge, and its mo- 
ments, are measurable and testable already at the LHC. 
With potential to uniquely determine quantum number 
of certain new physics particles, should they show up, it 
is important to verify jet charge on standard model pro- 
cesses. Thus jet charge holds promise as a measurable, 
calculable and useful observable. 
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